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Abstract 

This paper treats boundary conditions on black hole horizons for the full 
3+1D Einstein equations. Following a number of authors, the apparent hori- 
zon is employed as the inner boundary on a space slice. It is emphasized that a 
further condition is necessary for the system to be well posed; the "prescribed 
curvature conditions" are therefore proposed to complete the coordinate con- 
ditions at the black hole. These conditions lead to a system of two 2D elliptic 
differential equations on the inner boundary surface, which coexist nicely to 
the 3D equation for maximal slicing (or related slicing conditions). The over- 
all 2D/3D system is argued to be well posed and globally well behaved. The 
importance of "boundary conditions without boundary values" is emphasized. 
This paper is the first of a series. 
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I. INTRODUCTION 



If we need to define a flat surface, an elastic membrane — such as a soap bubble - 
stretched over a fixed rigid ring does a good job. The shape of the membrane is governed 
by the Laplace equation, with Dirichlet boundary conditions at its edge determined by the 
ring. What if we additionally need the membrane to have a hole somewhere near its center? 
If there is a small rigid ring available, we can simply drop it onto the membrane, and then 
cut out a hole in the membrane, fixing the new, inner edge to the small ring, again with 
Dirichlet boundary conditions. Overall, the small ring floats according to forces from the 
membrane, but the ring determines the size and shape of the hole. What do we do if no 
suitable small rigid ring is available? Our metaphor, admittedly loose, of course pertains 
to the study of black holes by numerical relativity, where indeed there is no suitable rigid 
"ring" . 

If we have closed loop of elastic string, it can serve as the inner boundary. The size 
and shape of the hole is not fixed, but rather determined by force balance between the 
string and the membrane. We may guess the equilibrium shape of the hole to be a circle. 
The curvature of the circle, and hence the size of the hole, will be determined by balance 
between tension along the string and surface tension in the membrane. The equilibrium will 
be described by some differential equations: In the membrane, the 2D Laplace equation; and 
along the string, two stationary wave equations — two, because the string has two transverse 
degrees of freedom in space.Q Boundary conditions at the edge of the hole couple all of these 
equations together. 

Let us now pass from the metaphorical soap bubble to the real black hole. Relativists 
have long hoped to pose numerical boundary conditions at the horizon of a black hole, 
and use the horizon as the inner boundary of the numerical grid. Recently, progress has 
been achieved in spherical symmetry by Scheel, Shapiro & Teukolsky , and by Anninos, 
Daues, Masso, Seidel, and Suen || (see also [|J). These two groups carried out important 
demonstrations of the feasibility of such horizon boundary conditions by evolving spherically 
symmetric black holes (in Brans-Dicke theory and Einstein theory respectively). The basic 
idea is to make the inner boundary an apparent horizon. 

It would be highly desirable to extend such horizon boundary conditions to more complex 
situations involving distorted or rotating black holes, and black hole binaries. The main 
purpose of this paper is to work out the underlying theory of these coordinate conditions 
in 3D generality, as a step of such an extension. In particular the 2D differential equations 
governing such inner boundary surfaces will be studied with a view toward making them 
well-posed, when connected to the 3D equation for maximal slicing in relativity, or to various 
other slicing conditions (i.e., choice of lapse) such as generalizations of maximal slicing [Q. 
The basic idea of using an apparent horizon will turn out to work, but several difficult 
issues crop up along the way, and must be dealt with. We need boundary conditions at the 
horizon, but we do not know how to supply boundary values there; how to resolve this issue 
of "boundary conditions without boundary values" is explained along the way. 



The real string also has a longitudinal degree of freedom; but this can be decoupled by attaching 
the membrane to the string with a slip boundary — as with a soap bubble. 
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The development in this paper will not make specific assumptions about the conditions 
for 3+ ID slicing and 3D spatial coordinates, though we will often assume that these condi- 
tions are implemented by some kind of 3D elliptic or parabolic equations, which will require 
boundary conditions at the hole. Hyperbolic operators (e.g., arising in harmonic slicing or 
harmonic coordinates) form quite a different case, and probably a simpler one, since they 
do not require boundary conditions at the hole.f] (I am grateful to Greg Cook for pointing 
this out.) 

In future papers in this series, further topics will be treated: 

• Horizon boundary conditions for stationary black holes, relevant to the late stages of 
numerical calculations. 

• The connection between spatial coordinate conditions, i.e., choice of shift, and horizon 
boundary conditions; this is especially important for rotating black holes. 

• Boundary conditions for the constraint equations of general relativity on horizons; 
again a key issue is "boundary conditions without boundary values" . 

II. MOST APPARENT HORIZONS ARE WILD 

Hawking || defined trapped surfaces and apparent horizon as follows. 

Definition 1. A trapped surface S is an achronal 2-surface in spacetime for which the 
outward null convergence obeys 

p>0 (1) 
A marginally trapped surface S is one for which the outward null convergence obeys 

p = (2) 

The basic property of these surface is: 

Proposition 2. (Hawking & Ellis fj/J. Given a trapped or marginally trapped surface S. 
No point p on S can lie outside the event horizon. 

Definition 3. In a given space slice V, an apparent horizon is the outward boundary of 
trapped surfaces that lie in V. 

Notice that the given slice V is an essential ingredient of the definition. To say that "We 
are going to choose a slice that meets an apparent horizon..." is a circular definition. This 
issue is not just a niceity; to ignore it can lead to serious trouble, as we will soon see. One 
way to see the problem is to observe that any slice passing far enough into the black hole - 
over a wide range — will meet an apparent horizon, and therefore will obey the boundary 
condition. A condition that excludes nothing is not a useful condition. 



2 On the other hand, it is not clear whether harmonic slicing avoids coordinate singularities. 
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To de-circularize the definitions, we must drop the slice, and refer not to an apparent 
horizon, but to a marginally trapped surface. So now we can sensibly say "We are going 
to choose a slice that meets a marginally trapped surface." And we can proceed to build 
a 3+1D code that uses marginally trapped surfaces S as the inner boundary of the slices. 
However, this code may not work; it is vulnerable to crashing after a short time. The pitfall 
is that "most" marginally trapped surfaces are wild surfaces; this is very unlike the familiar 
situation with maximal slices, which are automatically smooth (technically, thanks to elliptic 
regularity). The wildness of marginally trapped surfaces is intrinsically a non-spherically- 
symmetric phenomenon, and does not show up at all in spherically symmetric setting - 
consistent with the success of Anninos, Daues, Masso, Seidel, and Suen. 

A related property of a generic marginally trapped surface, is that it can always be 
deformed spatially outward from any point on it: 

Proposition 4. Given any smooth marginally trapped surface S, such that either a or 
T^ u l^l v does not vanish identically on S. Given any point p on S, and given a spacelike 
outward-pointing vector u at p. Then S can be locally perturbed into a 1-parameter family 
of marginally trapped surfaces S(e) so that p moves in the u-direction. 

(Here a is the outward null shear, T ilv lH v is the outward null component of the stress 
energy tensor, and we will assume throughout that obeys the dominant energy 

condition. Actually Proposition 4 is true for any u, but is most interesting if u is spacelike 
outward-pointing. ) 

This leads to a puzzle. Start with an apparent horizon S and some point p on it. What 
is to keep us from continually perturbing it, extending the 1-parameter family 5(e) , until 
p(e) passes outside the event horizon, contradicting Proposition 2? The only thing that can 
go wrong is that S must "go wild" — i.e., cease being smooth — first. 

In fact, the following conjecture, at first a bit startling, can be surmised as the obvious 
general answer to the puzzle. 

Conjecture A. The outward boundary in spacetime of marginally trapped surfaces is the 
event horizon. 

I have proven this conjecture under some assumptions, plausible but not rigorously es- 
tablished, about how spacetime settles down to a nonextremal black hole subsequent to 
gravitational collapse. The trapped surfaces S that pass close to the event horizon are wild. 
Let us give an example to illustrate what happens. 

Sometime during the dynamical phase of collapse, after the black hole has formed but 
before it has settled down, choose a slice that contains an apparent horizon S, and choose 
a point p on S. Also choose any point q on the intersection of the slice with event horizon. 
In the usual picture, q is considerably outside the apparent horizon, so it is unclear how 
any trapped surface can pass close to q.) However, by Proposition 4, we can now extend a 
1-parameter family of marginally trapped surfaces S(e), so that p(e) remains in the slice V 
(though most of S(e) does not) and so that p(e) — > q at some parameter value cq. What 
happens is that all of S(e) approaches the event horizon as e ^ e , but that most of S(e) 
moves far to the future. Therefore the trapped surfaces S which pass close to q look as 
follows. Most of S lies far to the future, very close to the event horizon of the settled-down 
black hole. Only a thin tendril of S extends back near q, in a thin tubular neighborhood 
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of a generator of the event horizon. Though this tendril is nearly null, it is still part of a 
spacelike surface. 

Let us see what this means for numerical relativity. Each marginally trapped surface 
S(e) can serve as the inner boundary of a slice V(e), say a maximal slice, extending to spatial 
infinity. Moreover we can arrange for e to be proper time at infinity, so that the slicing "goes 
wild" at finite time. This is the serious trouble that such a code can run into. 

Details will be published elsewhere. 



III. THE PRESCRIBED CONVERGENCE CONDITIONS FOR 2-SURFACES 

This is not to say that apparent horizons are necessarily a bad idea for inner boundary 
conditions, only that an ingredient is missing. We propose that the missing ingredient is 
the simplest possible thing, the convergence p' of the inward null normal to a surface S, and 
that the well posed way to put an inner boundary on a slice is to prescribe both the outward 
null convergence p and the inward null convergence p 1 . 



A. The Two Mean Extrinsic Curvatures of a 2-Surface in Spacetime 

Consider a spatial 2-surface S immersed in spacetime. (Throughout, S will be topolog- 
ically a 2-sphere unless otherwise noted.) All over S, we can choose an orthonormal frame 
of reference so that the time axis eg and one of the spatial axes ej are normal to S, 
while the other two spatial cLXGS €?2 ; fig £1X6 tangent to S. Then S has defined on it two mean 
extrinsic curvatures, namely H in the 6-direction, and Hi in the 1-direction. If ^/h is the 
element of area on S, then in suitable local coordinates we can take as the definitions 

d^Vh = —H^y/h (3a) 
diVh = -HiVh (3b) 

(the minus sign in these equations is a matter of convention). In any coordinate system, 
these mean extrinsic curvatures form a 4-vector H^ orthogonal to S. 

Alternatively we can follow null methodology |7J] and choose two orthonormal null vectors, 
the outward and inward null normals to S respectively, as 

^ eA vr (4b) 

with = = n^n^, 1^ = —1, and then the null convergences of 1 and n are 

(5a, 

' - ^ (5b) 
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1 


-> 1 = 


e^l 


n — 


-> n = 


e^n 


P ~ 


P = 


e^p 


P'~ 


>P~' = 





Notice, however, that we had to choose some frame of reference on S to define ifg and 
(or p and p'). The arbitrariness in choice of frame amounts to a boost in the 1-direction 
(a GHP boost) all over S 0. Under such a boost these quantities transform like 

(6a) 
(6b) 
(6c) 
(6d) 

where the boost rapidity if) is an arbitrary function on S. Only one boost-invariant scalar 
can be formed from the mean extrinsic curvatures, namely H^H^ = —H? + H? = —2pp'. 

How then do we determine a 2-surface S in spacetime? Since S has two transverse 
degrees of freedom, we need two conditions at each point of S, that is, two equations in two 
coordinates. It is tempting therefore to put conditions on the two quantities {Hq,Hi} or 
equivalently {p, p'}\ 

p = (some fixed function on S) (7a) 
p' = (some other fixed function on S) (7b) 

(For instance, to use an apparent horizon we can just take p = as one of the two conditions.) 
These equations will be called the equations of prescribed convergence, or, for short, the PC 
equations. However the PC equations are not enough by themselves, because of the boost 
arbitrariness, Eq. (|J). We are missing one condition, namely something to fix if> as a function 
of two coordinates on S. 

If we have already decided on a slicing condition for spacetime, e.g., maximal slicing, or 
something related to it, then the slice itself provides the missing condition: We can use the 
unit timelike normal to the slice to fix the O-direction at S. Then the PC equations, together 
with the maximal slicing equation, are just enough. This is of course no guarantee that the 
equations form a well-posed system, meaning that V and S can never "go wild" . However it 
is a reasonable conjecture that this is so, and this paper will present a considerable amount 
of evidence in favor of this conjecture. 



B. The 3+1 and 2+1 Splits 

The spacetime metric is 

ds 2 = -(a 2 - frp^dt 2 + 2f3 i dx i dt + % j dx i dx j (8) 

The notation is now necessarily going to become a little complicated, so the reader is 
asked to be patient. Starting from the full spacetime geometry, we choose slices and carry 
out the standard 3+1 split. The unit future-pointing timelike normal to the slices is denoted 
eg. Spacetime indices run over p, u, . . . = 0,1, 2, 3, where denotes time i; spatial indices 
run over = 1, 2, 3, and can be lowered and raised with the spatial metric 7^ and its 

inverse 7 lJ . 

Each space slice t = const has an inner boundary on a 2-surface S. For convenience in 
this paper, we will throughout choose spatial coordinates {x 1 } and shift vector (3 l so that 
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S always lies at r = tq = const, where r = x ; (9) 

This means that /3j is not freely specifiable at S, but must be chosen to match the motion 
of S: 

B = Pie\ fixed at S by the PC conditions (10) 

here e~ is the unit outward spatial normal to S. One would actually like more general spatial 
coordinates, to allow the black hole to move through the spatial coordinate system — to 
"fly through the grid" . The generalization to such coordinates is straightforward but is not 
developed in this paper. 

The Riemannian 2-metric on S will be denoted h ab ; surface indices run over a,b, . . . = 2, 3 
and can be lowered and raised with h ab and its inverse h ab . We use a 2+1 split at S, and 
2-tensors on S will carry indices a,b, . . .. The spatial extrinsic 2-curvature of S with respect 
to the spatial normal direction el is a 2-tensor 

H a b = —-£-ih a b (11) 
where C denotes Lie derivative, and the spatial mean extrinsic 2-curvature i7j is its 2-trace, 

Hi = h ab H ab = -^d^ (12) 

where h = deth a b- The Gauss-Codazzi equations of the 2+1 split imply 

3 R = 2d i H i - Hi - H ab H ab + 2 R (13a) 
h ab 3 R ab = dtHt - H ab H ab + 2 R (13b) 
2 D a K{ = diH 6 - H^ + H^yK (13c) 

where 2 R a b is the Ricci 2-tensor of h ab and 2 R = h ab R ab . 

We also need the timelike extrinsic 2-curvature of S with respect to the timelike normal 
direction eg; it is the 2-tensor gotten by 2+1 projection of Kij, 

Jab =± K ah (14) 

and the timelike mean extrinsic 2-curvature Hq is its 2-trace, 

#o = h ab J ab = h ab K ab ; (15) 
here + denotes the projection of 3-tensors into 2-tensors at S. 



C. Evolution of the Mean Extrinsic Curvatures 

We will present a system of equations which is fully covariant under transformations of 
the spatial coordinates x % and shift f3 l that preserve the inner-boundary constraints at S, 
Eqs. (|^, |i~0p . That is, we shall derive a system of equations which can be used with any 
choice of spatial gauge away from S. Some intermediate calculations, though, are most 
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easily done in a particular spatial coordinate system, namely, Gaussian normal coordinates 
in space. Thus our spatial metric can be taken as 



r y i jdx l dx J = dr l + h ab dx a dx b 



(16) 



where r = x l . The inner boundary in spacetime is the hypersurface r = r$ = const, 
consisting of all the S for all t, which we want to be a spacelike or null hypersurface lying 
within or on the event horizon. Values at S will be denoted " \ s " . Thus by Eq. (|T0D 



A I 



B 



and also we will denote the inner boundary value the lapse at S as 

a\ s = A 

and its spatial normal derivative as 

dials = M 



(17a) 
(17b) 
(17c) 



At the inner boundary surface S, the two mean curvatures H§\ s and Hi\ s are the quan- 
tities we wish to prescribe. Then we wish to derive these inner boundary values for the lapse 
a and shift f3 l . 

From these definitions and the Einstein equations follow the evolution equations for Hq 
and Hf 



1 -{H 2 - H 2 ) + \{H ab H ab + J ab J ab ) - K a r K r a + 8vrT 5Q + \ 2 R - 2 A 



+ 
+ 



J ab H ab + - H^iK - 8vrT 5r + ( 2 D a K?) + 2K a r 2 D a 



B 



Hi 



J ab H ab - 8vrT ar - ( 2 D a K?) - 2K a r 2 D a 



A 



+ 



-(H? + Hf) + ^(H ab H ab + J ab J ab ) + K a r K r a - H 6 tiK 



8nT< 



00 



- 2 R+ 2 A 
2 



B 



+ 



A 



(18a) 



(18b) 



where d[ n ^ denotes the projection of the time evolution operator d t = d/dt normal to S, 



«9 (n) = Ae 11 -^— + Be^-^— 
1 dxv 1 dx^ 



(19) 



These equations are valid in any 3+1 coordinate system that obeys Eqs. (|9|, [T(]) at S. 

If we now view Hq and Hi as prescribed functions of (t, x a ), Eqs. ( |18|) become a set of two 
2nd-order coupled partial differential equations in {x a } on 5* for three unknown functions A, 
B and Aj. Thus, the inner boundary conditions for the lapse a and the radial component 
of the shift /3 1 are fixed. The tangential boundary value of the shift is still free, and can be 
chosen separately to enforce, say, "slip" or "no-slip" conditions as desired. 
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D. The Naturality of the Inner Boundary Conditions 

It may seem surprising that these equations involve Ai as well as A and B\ however a 
little reflection shows this to be natural and in fact desirable. To illustrate this, restrict to 
maximal slicing trK = temporarily. The maximal slicing equation 

= ( 3 A - KijK^ - matter)a (20) 

is a 2nd-order elliptic equation in (x l ) for a, and therefore admits at the inner boundary S 
any of these well-posed boundary conditions: 

a\ s = A Dirichlet boundary conditions, or 
®i a \s = A-i Neumann boundary conditions, or more generally, 
Fdq>\ s + -F/v<9ja| 5 = homogeneous mixed boundary conditions; 
F D a\ s + F N dia\ s = F inhomogeneous mixed boundary conditions; 

where A, Aj, Fd, Fn, and F are prescribed functions of (t, x a ). 

Conditions (pT|cj21|d) simply say that there exists a prescribed linear relation between a 
and aj on S, but that neither is fixed individually; the distinction between these conditions 
(pl|c,pl|d) is the function F on the right-hand-side, which may be either (homogeneous 
case) or a given function (inhomogeneous case). 

Now, given Hq and iff, Eqs. ( |T8| ) do give two prescribed linear relations among the three 
functions (B — A), (B + A) and A^, albeit implicit ones. One can imagine using one of 
Eqs. flTE| ) to eliminate B, whereupon the other becomes a single linear relation, in the form 
of Eq. (pl|c), between A and A^ — albeit an implicit one, involving Green functions of the 
operators V and V . Thus, Eqs. ( |i"8| ) do appear to give well-posed boundary conditions 
for maximal slicing, Eq. (|20"D, and they also should do so for attractive generalizations of 
maximal slicing ||. 

In fact, not only is it admissible to have homogeneous mixed boundary conditions for 
the lapse a at S, it is also desirable. Imagine we instead used Dirichlet boundary conditions 
(|T7|a). How big do we then make the inner boundary value AI — by which we mean, what 
fixes the overall scale of A at the inner boundary? We would like to make it big enough so 
that, over a long evolution, the inner boundary just "keeps up" with the outer boundary, 
neither shooting way ahead, or falling way behind. How big is that? The first guess is 
A — 1, but that cannot be right, because the inner boundary is in the strong field region. 
The lapse equation is not going to tell us how big to make A, precisely because it is happy 
with any A. Thus a pure Dirichlet boundary condition is not generally going to work well 
for long evolutions. 

Thus, mixed homogeneous boundary conditions resolve this issue. They provide "bound- 
ary conditions without boundary values" . 

A homogeneous mixed boundary condition can work well — it can be thought of as a 
"feedback mechanism" . If the inner boundary falls behind the outer boundary, then A^ will 
become large and positive, and the mixed boundary condition then can "tell" A to become 
larger. If the inner boundary shoots ahead of the outer boundary, then Ai will become large 
and negative, and the mixed boundary condition then can "tell" A to become smaller. A 
concrete example of this "mixed-boundary-condition feedback mechanism" will be presented 
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(21a) 
(21b) 
(21c) 
(21d) 



below. Since the boundary condition is homogeneous, no data need be given at the inner 
boundary to fix the overall scale of A; this scale is self-adjusting. 



E. The PC Equations — Inner Boundary Conditions for Prescribed p and p' 

From Eqs. (§,0) can be found the evolution equations for the null convergences p and 
p'. When we take p and p' as prescribed functions of (t,x a ), these equations become the 
equations of prescribed convergence, or, for short, the PC equations: 



d(n) 
t P 



1 

72 



p(p' + 2p- trK/y/2) + -£ 



+ 



dt P 



1 

V2 

P 
1 " 

i 

p' 



p(-p' + p- trK/y/2) + -V 

A 

f/(- p + p> - trK/y/2) + \v' 
p\p + 2p' - trK/V2) + -£' 

M 



(B + A) 
IB -A) 



(B + A) 
IB -A) 



(22a) 



(22b) 



a set of two 2D elliptic equations relating three unknown functions A, B, Aj, where 



V = 


h ab CD a + K r a )CD b + K r b ) 


2 






v' = 


h ab CD a -K r a )CD b -K r b ) 


2 






£ = 


\(Jab-H ab )(J ab -H ab ) 


-2p 2 + 




-Tor) 


£' = 


\{J ab + H ab ){J ab + H ab ) 


- 2p' 2 + 


8tt(T 58 


+ T 0r ) 



(23a) 
(23b) 
(23c) 
(23d) 



The latter quantities obey some useful relations. The two differential operators T> and T>' 
live on S; they are not generally self-adjoint, due to the terms in K r a1 but are the adjoints 
of each other: 



= D'. 

The two scalar functions £ and £' are nonnegative: 

£ > 
£' > 



(24) 



(25a) 
(25b) 



as long as matter obeys the dominant energy condition, which we will assume throughout. 
These equations are valid in any 3+1 coordinate system that obeys Eqs. (|||, [1^) at 5*. 
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F. Special Case: Maximal Slices Bounded by Apparent Horizons 



As a special but important choice of coordinate conditions — still general enough to 



study binary black hole coalescences in 3+1D — let us take: 

trK = (maximal slicing) (26a) 

p = (apparent-horizon inner boundary) (26b) 

p = F(x a ) (some prescribed function) (26c) 

Then the PC equations reduce to 

= £{B + A) +V(B - A) (27a) 

= {V + 2p' 2 )(B + A) + (£' + Ap' 2 )(B - A) + 2v / 2p'A i (27b) 



G. Special Case: Stationary Black Holes 

Stationary black holes are relevant because calculations of binary coalescence will eventu- 
ally settle down to a stationary black hole, and coordinate conditions are desirable that will 
somehow "lock onto" the geometry of the stationary black hole and render it recognizable. 
For a stationary black hole, we have 



and then from Eq. 



p = 
8 = 

) we have immediately on S 
B = A 



(28a) 
(28b) 



(29) 



so that the PC conditions reduce to one 2D elliptic equation relating two unknown functions 



A, A r . 



d[ n) p' = -V2 



p\p' - trK/V2) + -V 



A - 



P 



(30a) 



The operator on A is not obviously invertible; however in a later paper in this series, it 
will be proven that this operator is in fact invertable under fairly general conditions for a 
nonmaximal stationary black hole, and so the PC conditions can be expected to successfully 
"lock onto" the Kerr black hole at late times. 



IV. CONCLUSION 

We have studied the apparent horizon boundary condition p = on a spacelike 2-surface 
S as an inner coordinate condition at a black hole. Since most apparent horizons are wild, 
another condition is required to ensure a well posed evolution. We propose the Prescribed 
Curvature Equations, or for short the PC Equations 
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p = (some fixed function on S) 
p' = (some other fixed function on S) 



(31a) 
(31b) 



an effective condition at the black hole. These equations have the following properties: 

• They give rise to a system of two linear 2D elliptic equations on S for three unknown 
boundary values of the lapse and shift. 

• They therefore connect well with maximal slicing, and with related slicing conditions 
that involve 3D elliptic or parabolic equations. 

• They are well posed and solvable for stationary black holes. 

• They appear likely to be well posed and solvable under very general conditions. 

Some additional numerical work will be required to solve the PC equations as part of a 
numerical relativity code. However, solving these equations cannot be not much harder than 
finding apparent horizons in the first place, (see, e.g., and may well be easier because 

these equations are linear. Therefore the extra work seems unlikely to be prohibitive. 
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